AD-R16e  976  FIXED  POINT  HETHOOS  FOP  THE  COHPLENENTRRITV  PROBLEM(U)  1/1 
UISCONSIN  UNIV'HAOISON  HATHEMOTICS  RESEARCH  CENTER 
P  K  SUBRANANIAN  AUG  8S  HRC-TSR-2857  DAAG29-88-C-8e41 

F/6  12/1  NL 


UNCLASSIFIED 


K'- 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 
610  Wainut  Street 
Madison,  Wisconsin  53705 


August  1985 


(Received  August  5,  1985) 


1C 


i 

ECECTE 
NOV  7  1985 


Approved  for  public  release 
Distribution  unlimited 


Sponsored  by 


yCj 


U.S.  Amy  Research  Office 
P.O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


National  Science  Foundation 
Vlashington,  D.C.  20550 


85  11  06  0 


- 

UNIVERSITY  OF  WISCONSIN-MADISON 
MATHEMATICS  RESEARCH  CENTER 

FIXED  POINT  METHODS  FOR 
THE  COMPLEMENTARITY  PROBLEM 


P.  K.  Subramanian 


Technical  Summary  Report  2857 
August  1985 

ABSTRACT 

This  paper  is  concerned  with  iterative  procedures  for  the  monotone  com- 
plementarity  problem.  Ow  iterative  methods  consist  of  finding  fixed  points 
of  appropriate  continuous  maps.  In  the  case  of  the  linear  complementarity 
problem,  it  is  shown  that  the  problem  is  solvable  if  and  only  if  the  sequence 
of  iterates  is  bounded  in  which  case  summability  methods  are  used  to  find  a 
solution  of  the  problem.  This  procedure  is  then  used  to  find  a  solution  of  the 
nonlinear  compler  lentarity  problem  satisfying  certain  regularity  conditions 
for  which  the  problem  has  a  nonempty  bounded  solution  set. 

AMS(MOS)  Clkssification:  90C30,  90C25 
Keywords:  Monotone  operators,  Complementarity  problem 
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dation  Grants  DCR-8420963  and  MCS-8102684. 


SIGNIFICANCE  AND  EXPLANATION 


Fixed  point  iterative  procedures  are  defined  which  use  summability  the¬ 
ory  for  the  solution  of  the  monotone  complementarity  problem.  Mangasarian 
^  McLinden  have  established  the  existence  and  boundedness  of  solution  set 
under  certain  regularity  conditions.  In  this  paper  we  provide  a  constructive 
procedure  to  obtain  a  solution  under  these  conditions. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC.  and  not  with  the  author  of  this  report. 


FIXED  POINT  METHODS 


FOR  THE  COMPLEMENTARITY  PROBLEM 
P.  K.  Subramanian 

1.  Introduction 

We  are  concerned  in  this  paper  with  the  complementarity  problem,  viz.,  that 
of  finding  a  Zo  (if  it  exists)  such  that  F{  Zo)>0  and  such  that  zj F(zo)  =  0. 
Here  F  is  an  operator  from  9?”  to  3?”  .  In  particular,  we  are  concerned 
with  the  case  when  F  is  monotone,  that  is 

(i  -  y)^(F(x)  -  F(y))  >  0  Vx,ye»”. 

The  operator  F  is  strongly  monotone  if  there  exists  a  positive  real  number  A 
such  that 

(x  -  y)^(F(x)  -  F{y))  >  Ajjx  -  yf . 

When  F  is  an  affine  map,  F(x)  =  Mx+q,  we  shall  refer  to  the  complementar¬ 
ity  problem  as  the  linear  complementarity  problem  and  write  LCP{M,  q)  in 
this  case.  Otherwise  we  shall  refer  to  the  complementarity  problem  as  the 
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In  the  case  of  LCP{M,  q) ,  when  M  is  positive  semidefinite,  if  the  prob¬ 
lem  is  feasible,  that  is  there  exists  a:  >  0  such  that  Mx  +  9  >  0,  the  problem 
is  solvable  [Eaves,  1971].  This  is  not  the  case  for  NLCP{F)  ([Megiddo, 
1977],  [Garcia,  1977]).  However,  for  e  >  0,  if  we  consider  the  Tihonov  reg¬ 
ularization  Fc  \=  F  +  el  then  the  corresponding  problem  N LC P(Fe)  has  a 
unique  solution  since  F*  is  strongly  monotone  [Karamardian,  1972].  When 
e  — >  0,  then  x*,  the  solution  of  NLCP{Fe),  converges  to  the  least  two-norm 
solution  of  NLCP(F)  ,  provided  NLCP{F)  is  solvable  [Brezis,  1973]. 

A  solution  of  NLCP{F)  is  also  a  fixed  point  of  the  map 

I  «  (z-F(x))^  :=  max  {0, x  -  F(x)}, 

The  principal  aim  of  this  paper  is  to  consider  iterative  procedures  to  find  such 
fixed  points.  We  shall  show  that  in  the  linear  case  the  sequence  of  iterates 
is  bounded  if  and  only  if  LCP{M,  q)  is  solvable.  When  this  is  the  case,  we 
use  summability  methods  to  obtain  a  solution  of  the  problem.  Although 
feasibility  of  the  monotone  NLCP{F)  does  not  imply  its  solvability,  it  is  a 
theorem  of  Mangasarian  ic  McLinden  [1985]  that  when  a  regularity  condition 
such  as  the  distribute  Slater  constraint  qualification  is  satisfied  then  the 
solution  set  is  bounded.  We  show  how  the  iterative  procedure  for  the  linear 
case  may  be  adapted  to  find  a  solution  in  this  special  case. 

We  briefly  describe  the  notation  used  in  this  paper.  We  use  for  the 
space  of  real  ordered  n-tuples  .  All  vectors  are  column  vectors  and  we  use  the 
Euclidean  norm  throughout.  Given  a  vector  x,  we  denote  its  component 
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by  Xi.  We  say  i  >  0  if  Xi  >  OVi.  The  nonnegative  orthant  is  denoted  by 

We  use  superscripts  to  distinguish  between  vectors,  e.g.,  x',z^  etc.  For 
X.  y  *?  3?"  x^  indicates  the  transpose  of  x,  x^y  their  inner  product.  Occa¬ 
sionally.  the  superscript  T  will  be  suppressed.  All  matrices  are  indicated  by 
upper  case  letters  A,  B,  C  etc.  The  row  of  A  is  denoted  by  A,  while  its 
y'^column  is  denoted  by  A.j.  The  transpose  of  A  is  denoted  by  A^. 

Given  NLCP{F)  .  we  define  the  feasible  set  and  solution  set  by  5(F) 
and  5(F)  respectively,  that  is, 

5(F)  =  {x£  :  F(x)  e  ^Stl) 

5(F)  =  {x  6  5(F)  :  x^F(x)  =  0}. 

In  the  case  of  LC P(M,  y ) ,  we  shall  denote  these  sets  by  5 (A/,  q)  and  5(A/,  q) 
respectively.  Finally  the  end  of  a  proof  is  signified  by  |. 

2.  Fixed  point  methods 

We  begin  with  the  well  known  notion  of  a  contraction  mapping. 

2.1  Definition.  Let  P  :  D  C  5R”  — >  3?”  .  We  say  P  is  Lipschitzian  with 
modulus  L  >  0  if 

\'P{^)  -  ^(y)li  <  L\\x  -  yl!  Vx,  y  €  D. 

When  L  <  \  {L  <  1)  «>e  say  P  is  non-expansive  (contractive)  . 

The  following  Theorem  is  classical;  see  e.g.,  [Ortega  and  Rheinboldt, 
1970,  page  120]. 
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2.2  Theorem.  (Banach’s  contraction  mapping  princplt).  Let  P:  D  C 
3?"  — ►  3?”  and  let  Do  be  a  closed  subset  of  D  such  that  PDo  =  {/’(at)  : 
X  €  Do}  C  Dq.  If  P  is  a  contraction  mapping  on  Do  with  modulus  L,  then 
P  has  a  unique  fixed  point  x  in  Do-  Further,  for  any  point  x^  in  Do,  the 
sequence  {x*}  where  =  P{x*),  converges  to  x  with  the  following  linear 
rate  : 


The  content  of  the  following  proposition  is  well  known.  We  state  it  in 
following  form  for  later  use  and  furnish  a  proof  for  the  sake  of  completeness. 
2.3  Proposition.  Let  F:  D  C  3?"  — ►  3?"  be  monotone  and  Lipschitzian 
with  modulus  L.  Suppose  that  e  >  0,  a  >  0  and  ea  <  1.  Then  the 
projection  map  IP  defined  by 

lP(x)  =  |x  -  a(F(x)  +  ex)  |  ,  x  €  D 

is  also  Lipschitzian  with  modulus  k{a)  ~  v^(l  -  ae)^  -|-  (oL)*.  If  a  < 
2e/ y/s^  +  L^,  then  IP  is  a  contraction  and  k  attains  its  minimum  value 

fcmin  (o)  =  L/  \/L2  +  e2  for  Q  =  e/L^  +  e^. 

Proof 

We  have 

l!IP(x)  -  lP(y)||^  =  ll{x  -  a(F(x)  +  ex)}^  -  {y  -  o(F(y)  +  ey)}^||^ 

<  ||{x  -  a(F(x)  +  ex)}  -  {y  -  a(F(y)  +  ey)}f 


since  projection  on  9?"  is  non-expansive.  Hence, 


l|]P(x)  -  IP(y)!i^  <  ||(x  -  y)(l  -  eo)  -  a(F(i)  -  F{y))\\^ 

=  ilx  -  y||^  (1  -  ae)^  +  a^\  F{x)  -  F(y)!:^ 

-  2a(l  -  ae)(i  -  y)(F(i)  -  F(y)). 

Since  ae  <  1  and  <  F  '(:r)  -  F(y),  X  -  y  >>  0  from  the  monotonicity  of  F, 
|iIP(i)  -  IP(y)i|^  <  llx  -  yll^  {(1  -  a£r)2  +  {aF)^}. 

The  other  claims  about  k{a)  are  easy  to  verify.  | 

2.4  Theorem.  Let  F  :  IR”  — »  !R”  be  monotone  and  Lipschitzian  with 
modulus  L.  Let  {cn}  be  a  sequence  of  positive  reals,  i  0,  For  n  =  1,  2,  •  •  • 
let 

IP„(a:)  =  {l  -  Qn  (F(x)  +  £nx)  }  ^ 
and  for  m  =  1,  2,  •  •  •  and  x  €  3?"  let 


ip;['(x)  =  Pn  o  •  •  •  o  IPn  (j)  =  x(n,m). 

m  times 


Suppose  further  that 


For  n  =  1,  2,  •  •  • ,  let  x'^  be  defined  by 

i’‘ -  x(n,m),  where  ||x(n,m  +  l)  -  x(n,m)||  <  Sn- 


Then  the  sequence  {||x”||}  is  bounded  if  and  only  if  NLCP{F)  is  solvable 
and  in  this  case,  ^  x,  the  least  two-norm  solution  of  NLCP[F)  . 

Proof 

From  Proposition  2.3,  IP„  is  a  contraction  with  modulus  <  1.  By  the 
contraction  mapping  principle,  given  any  x". 

lim  IPn-'  — >  2",  P„(2”)  =  2". 

J  — »oo 

Note  that  2”  solves  NLCP{F  +  e„/)  uniquely.  Since,  by  definition, 

IP„(x(n,m))  =  i(n,  m  +  l) 


we  have 


6n  >  ||x(n,m  +  1)  —  x(n,m)||  >  ||x(n,m)  -  2”||  —  ||x(n,m  +  1)  -  2” 
and 


||x(n,m  +  1)  -  2"!!  =  j|IP„(x(n,m)  -  Pn(2")||  <  ■  |!i(n,m)  -  2”|i 

it  follows  that 

Sn  >  llx(n,m  +  1)  -  x(n,m)||  >  (1  -  A:„)  |jx(n, m)  -  2”|1 


and  that 
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The  conclusions  about  {x"}  follow  from  [Brezis,  1973]  (see  also  [Subrama- 
nian,  1985]).  | 

We  remark  that  the  last  Theorem  is  a  two-step  process  in  the  sense  that 
for  a  given  £„•  the  contraction  IP„  is  iterated  m  times  until  x(n,m)  is  close 
enough  to  the  solution  2"  of  NLCP(F  +  £n)-  One  then  takes  a  smaller 
£„  and  the  process  repeats.  Our  aim  now  is  to  prove  convergence  for  an 
algorithm  which  combines  both  steps  into  a  single  step.  We  shall  need  the 
following  notions  from  the  thet.ry  of  summability. 

2.5  Definition.  An  infinite  matrix  A  —  (Aij),  i,  j  =  1,  2,  •  >  ■ ,  is  said  to  be 
convergence  preserving  if  for  any  sequence  {in},  the  sequence  {t/n}  defined 

by 

00 

Vn  —  ^  ^  AnjXj  (2 -6) 

is  well  defined  and  limXn  =  limy„.  We  call  (j/n)  the  A-transform  of  {in} 
and  write  yn  —  y4({in}). 

The  following  Theorem  is  classical.  Its  proof  may  be  found  for  instance 
in  [Peyerimhoff,  1969]. 

2.7  Theorem.  (O.  Toeptlitz).  An  infinite  matrix  A  =  (/I,;  ),  i,  j  —  1,  2,  •  •  • 
is  convergence  preserving  if  and  only  if 

( 1 )  EJl  1  ■  ^ exists, 

{a,}  is  bounded, 

(2)  lim,(y;^~  1 /l,j)  =  1. 

(3)  lim,  Aij  0. 
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We  are  now  ready  to  prove  the  principal  theorem  of  this  paper. 

2.8  Theorem.  Let  M  be  a  positive  semidefinite  matrix.  Assume  that  that 
the  sequences  {on}  o,i^d  {cn}  of  positive  reals  are  such  that 

OO 

On  diverges, 


n  =  1 
OO 


converges, 

n  =  \ 

OO 

an^n  converges  and 


n=l 


(2.9) 


en  <  1,  I  0. 

Suppose  that  k  is  the  smallest  positive  integer  satisfying 


V^Pk  +  2pjt  <  L 


1 


(2.10). 


1  +  ||M|| 

Let  B  =  [Bij]  be  the  infinite  matrix  whose  n^^  row  Bn  is  defined  by 
D  _  f^l+k  02+k  an+k  ^  \ 

’°’"7 

where  Sn  =  =  i  Sy  +  *.  Let  =  0  and  having  x”,  determine  /rom 

^.n+l  _  _  an  +  k£n-\-k)x’^  -  an+k{Mx’^  +  9)}^  (2.11). 

Lc<  {y”}  be  the  B-transform  of  {x”},  that  is 


(2.12). 


Then  S(M,  g)  #  0 


{i”}  is  bounded.  When  this  condition  holds, 

yn  _  y*  ^  S(M,  g). 
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Proof 

Assume  that  k  satisfies  (2.10).  For  notational  convenience,  we  shall  write 

+  ^n+k  and  pri  Pn-rk‘ 

Obviuosly,  the  sequences  {^n})  {^n}  and  {Pn}  also  satisfy  the  conditions 
(2.13).  We  shall  write 


Fx  =  Mx  -i-  q,  FnX  =  Fx  +  SnX. 


Thus  we  can  write  (2.11)  in  the  form 

=  {(1  -  a„e-„)i"  -  OnFx"}^  (2.13). 

We  first  assume  that  {x”)  is  bounded  and  show  that  in  this  case  y”  — > 
y*  G  5(M,  q)  so  that  5(M,  q)  ^  0. 

Assume  then  that  {i”}  is  bounded.  Clearly.  3Ki  >  0  and  K2  >  0 
such  that 

llx"ll  <  Ki. 

||F„x"|l  =  ||Mx"  ^  q  +  e:nx"|| 

<  (1  +  IlMll)  •  |lx";l  ^  \\q\\ 

<  (1-  |;M!|)-A’,  +II9I1 

:=  K2. 
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Let  I  e  9?”  be  arbitrary  but  fixed.  Then  from  (2.13)  we  have 

<  ||(x”-x)-an(Fx"  +e„x")||2 

<  l|(x"-x)ii2-  2a„(Fx'^)(x"-x) 

-  2o„e„x”(x"  -  x)  +  (2.14) 

Since  M  is  positive  semidefinite  we  also  have 

(Fx’‘)(x"  -  x)  >  (i^x)(x''  -  x). 


Let 

A'3=sup  ljxl.(jx”-x||. 

n 

From  (2.14)  we  now  get 

2an(/’x)(x"  -  x)  <  ||x"  -  xf  -  ||z"+‘ -  il|2  +  2Q„£r„A3  + 

Summing  this  from  1  to  A:  we  obtain 

k  k 

2(Fx) an{x^  - x)  <  ||x‘-xf -||x^+‘-xf +  2A3  ^  a„Ern  +  A’| 

»J  =  I  n=l  n=l 

Divide  this  last  inequality  by  Sk  and  let  fc  — »  00.  From  the  assumed  prop¬ 
erties  in  (2.9)  of  the  sequences  {on},  {^n}  and  from  the  definition  of  {y"}, 
we  now  have 

lim  inf  <  Fx,  i  -  y*^  >  >  0. 

k 


•w'r:  w  ^ 


•  j»  •  V  •  \  *  •  •  •  *_>  *■  -  '  •  *  -  •  V  •  ji ^  *  I'  *1  Tji  W  *  _*•»_*  •  V.  •  fc  *  •  ■  ^  r  U*»  j  «  •%  W V  w  ■  w  J" 
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Since  y”  is  a  convex  combination  of  I^  it  follows  that  {i”} 

bounded  =>  {y"}  is  bounded.  Hence  {y"}  has  a  limit  point  y'  for  which 

<  Fx,  X  —  y*  >  >  0. 

Since  x  €  was  arbitrary,  y''  solves  LCP{M,  q) .  This  completes  our  proof 
that 

{x"}  bounded  ■=^  S{M,  q)  7^  0. 

Next  we  prove  y”  — >  y*. 

Since  S(M,  q)  7^  0,  choose  2  e  S{M,  q)  arbitrary  but  fixed.  By 
Theorem  2.4.2,  2  satisfies 

<  Fx",  x"  -  2  >  >  0  (2.15) 

since  x”  >  0.  From  (2.13)  and  (2.15)  we  have 

ll^n+i  _  2;||2  <  -2)-  a„(Fx”  +  e„x")f 

<  lix”  -  2|p  -  2an(Fx")(x"  -  2) 

-2ocnenx'^{x^-z)  +  cxlKl 

<  ||x"  -  2l|^  +  2Q„en  lx”(x”  -  z)\  -r  «n^2- 

Define  Pn{z)  by 

Mz)  :=  2c,„£„  ||i"||  III"  -  ^11  +  a;  (2.16) 

and  we  now  have 

-2||2  <  |jx--2||2+^„(2)  (2.17). 


mmm 

m 

HSE 

m 


Let  S  denote  5(M,  q)  and  let 


=  PeCx"). 


We  are  going  to  show  that  Bz*  such  that 


2*,  y" 


From  (2.17)  and  the  definition  of  e". 


n+l_,n+l||2  <  ||3;«+>  _  2"jj2  <  ||x"  -  + /?„(«). 


I"  '  -  Z 


Since  0n(^)  converges,  by  [Cheng,  1981,  Lemma  2.2.12],  we  can  conclude 


z'‘-2“||  converges. 


(2.18) 


By  parallelogram  law,  for  m  >  0, 

||^«+m  _  ^n||2  ^  2||x'‘+"‘  - 

2 

Since  5  is  convex,  (2”  +  2”''''”)/2  €  5.  Also,  2’*'^'”  is  the  closest  point  to 
x»i+rri  jjj  g  Hence, 


l^n+m  _  ^n||2  <  2||x"+'"  -  2'*||2  -  -2||x"+”*  -  2”+"*f . 


(2.19) 


Letting  z  =  z”^  in  (2.17)  and  noting  that  2"  is  the  closest  point  to  x"  in 
5,  it  follows  that  /3n(x’*)  <  /^n(x)  •  Now  let  2  =  2"  in  (2.17)  and  use 
induction  to  get 


m  >  0. 


Substitute  this  in  (2.19)  and  we  have 


|j^n  +  m  _  2n|j2  <  2\\x^  - 

n+  m 

-  4  2^  I3j{z).  (2.20) 

}=n 

From  (2.18)  and  the  fact  0n(z)  converges,  we  have  by  letting  n,m  -*  oo 
in  (2.20)  that 

ll^n+m  _  ^n||  ^  q 

SO  that  {2"}  is  Cauchy.  Since  S  is  closed,  3  2*  £  5  such  that  2”  —►2*. 
We  shall  now  show  that  y”  2’  as  well. 

Since  {y”}  is  also  bounded,  let  y*  be  any  of  its  limit  points.  Assume 
that  the  subsequence  y'**  converges  to  y*  .  From  our  proof  earlier,  y*  €  S  . 
Observe  that 


Notice  however  that 


is  simply  a  subsequence  of  the  B-transform  of  {z**}  ,  that  is  of 

{{”}  = 

Since  B  satisfies  all  the  conditions  of  Theorem  2.7,  it  is  a  convergence  pre¬ 
serving  matrix.  However,  2"  — >  z*  so  that  both  and  also  converge 
to  2*  .  If  we  take  limits  as  /c  — >  00  in  (2.22),  we  get 

<  y*  -  z*,  y*  -  z*  >  <  0 

so  that  y'  =  z* .  But  y*  was  any  arbitrary  limit  point  of  {y”} .  Hence 
y"  — »  2*  .  This  completes  our  proof  that 

{x”}  bounded  ==>  5(M,  9)  #  0  and  y"  — ►  z*  6  5(M,  g). 

We  now  prove  the  converse,  that  is  we  shall  assume  that  S{M,  9)  ^  0 
and  show  that  {x”}  is  bounded. 

Recall  from  (2.10)  that  k  satisfies 

+  <  L  =  j  ^  ||^j|- 

Hence  there  exists  a,  0  <  o  <  1/2  for  which 

yWk  +  Pk  <  (i+^)(,  +  |jM||)  ■ 


The  function  /(r) , 
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is  strictly  increasing  in  [0,  oo],  \imr  f{r)  -  L^r-  Thus  3f  >  0  such  that  for 
r  >  f, 

<  /(r)  <  /(r). 

Since  |  0  and  Pn+k  ~  have  for  all  n  >  0  and  r  >  f, 

V^  +  2p„  <  /(r)  <  /(r).  (2.23) 


By  assumption,  5^0.  Let  z  =  ij(0) ,  that  is  z  is  the  least  two-norm 
solution  of  LCP{M,  q) .  Define 

r  =  max(f,  -ik||)  +  1. 
o 

Our  aim  is  to  show  that 


-  zl!  <  r, 


Vn  >  0, 


that  is  {x"}  is  bounded  and  this  would  complete  our  proof. 

We  use  induction. 

For  n  =  0,  -  z\\  =  ||x||  <  ar  <  r. 

Suppose  now  that  ||z"  —  z\\  <  r .  Let  =  |jx"  -  z\\ .  From  (2.13), 


<  ilx"-x|i2-2a„e„x”(x"-x) 

-  2a„(Fx'')(x"  -  e)  -h  o^llFx”  enx'^W^ 


(2.24) 
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Since  z  e  S,  (fx'‘)(x'‘  -  z)  >0.  Also  if  Mn+i  <  Mn*  we  are  done.  So 
assume  /in+i  >  Mn  •  From  (2.24)  we  thus  get 

2an£nx^{x^  -  z)  <  a2||Fx'*  +  £r„x"ll2, 

that  is 

x^ix'^-z)  <  ^||Fx'‘  +  e„x'‘||2.  (2.25) 

Since 

l|x“||  <  III"  -  2||  +  ||z|| 

<  r  +  or 
=  (1  +  a)r, 

we  have 

||Fi"+£„i»||  <  ||Mi"+,||  +  £„||i"|| 

<  1IM||  ||i"||  +  ||«||  +  ||x»|| 

<  (t  +  <r)r-(l  +  l|M||)+l!«|| 

_  r 

~m 

=: C  say. 

From  (2.25)  we  now  get 

x"(i"-x)  < 

(X"-£)(X"-X)  <  -£(!»-£) 

<  y«"  +  W-|k"-*ll 


(2.26) 


Rewriting  this  last  inequality. 


Mn  <  ^t^^ronn. 


whence 


2fjil~2rofin  -  <  0. 


Since  f^n  >  0 ,  we  must  have 


P-n 


2ra  +  \/4r‘^o^  +  ^i^Pn 
4 

ro  (  2ro  +  2^y/2pn  \ 


and  finally 


Pn  <  r<T  +  |\/^- 


(2.27) 


Again  from  the  definition  of  in  (2.13), 


Pn+i  =  -  2| 


<  !|x'‘-z-a„(Fnx")ll 

<  +  a„||Fnx”|i 


^  Pn  Pn^^ 


(2.28) 


where  we  have  used  (2.26)  and  the  fact  q„  <  Pn  •  If  we  use  our  estimate 
of  Hn  from  (2.27)  in  (2.28)  we  get 


Pn+l  <  ra  ^\'2pn  +  Pnt 
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Substituting  for  ^  from  (2.26)  and  using  (2.23)  we  finally  get 

„  .  I  (v/2^  +  2pn)  r 

<  ra+  2 

.  .  /(O  ^ 


<  ra  + 


2  fir) 


since  a  <  1/2.  Hence  Hn+i  <  f  •  This  completes  our  induction  and  also 
the  proof  of  the  Theorem.  | 


2.29  Remark 

Our  proof  showing  that  {y"}  converges  by  considering  —  Pj(x")  is 
patterned  after  [Baillon,  1975],  who  uses  this  technique  to  construct  fixed 
points  of  non-expansive  maps.  Notice  also  Baillon ’s  use  of  the  Cesaro  matrix 
C  where  Cij  =  l/i  for  j  <  i.  while  Cij  =  0  for  j  >  i. 

3.  Application  to  NLCP(F) 

We  shall  now  show  that  the  proof  of  Theorem  2.8  can  be  used  to  construct 
a  solution  of  NLCP(F)  when  F  is  monotone  and  satisfies  some  regularity 
conditions  such  as  the  distributed  Slater  constraint  qualification  [Mangasar- 
ian  and  McLinden,  1985]. 

3.1  Definition.  Let  F:  D  C  — ►  IR"  .  We  say  that  F  satisfies  the 

distributed  Slater  constrainst  qualification  (DSCQ)  if  there  exist  p  points 
,■  ■  ■  ,,zf'  e  D,  nonnegative  weights  Xi,  X2,- ••  hp  (X^j  Ay  =  l)  such 
that  z  =  ^y  XjZ^  >  0  and  w  =  Ayu.’-^  >  0  where  ~ 


vv 


Mangasarian  and  McLinden  have  proved  the  following  Theorem. 

3.2  Theorem.  Let  F:  D  C  SR”  — SR”  ,  SR"  C  D  and  suppose  that  F  is 
monotone  and  continuous  on  D .  Assume  that  F  satisfies  (DSCQ).  Let 

p 

')  >  max  (l,  -zw  +  XjZ^w^), 

C  —  {z  C  SR"  :  wz  <  tbz  +  ')} 

where  Aj .  z  and  w  are  as  in  (DSCQ).  Then  N LC P{F)  is  solvable 

and  has  a  solution  z*  such  that  wz*  <  wz  + 

We  shall  now  show  that  the  technique  used  in  the  proof  of  Theorem  2.8 
can  be  used  to  construct  a  solution  of  N LC P{F)  guaranteed  by  Theorem 
3.2. 


3.3  Theorem.  Assume  that  F  satisfies  the  hypotheses  of  Theorem  S.2  and 
let  C  be  the  compact  convex  set  as  defined  in  that  Theorem.  Let  =  0 
and  given  i"  find  from  r"'*'* 


=  Pc{x^ 


n  ' 


Let  B  be  the  Cesaro  matrix  with 


and  let  {y"}  =  .  Then  y"  converges  to  a  solution  of  NLCP{F)  . 


Proof 


We  shall  give  only  a  brief  outline.  Since  {x"}  and  hence  {y"}  are  both 
bounded,  {y"}  has  a  limit  point  y* .  One  uses  the  monotonicity  of  F  to 
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show  that 

{F{y*,  x-y*)  >  0,  Vi  e  C. 

Hence  y*  is  a  fixed  point  of  the  map  x  •->  Pc{i  ~  ■^(^))  •  However,  Man- 
gasarian  and  McLinden  show  that  any  such  fixed  point  satisfies  wy*  < 
lyz  +  -7 .  Hence  y'  solves  NLCP{F)  .  One  can  now  show  that  y"  — >  y* 
by  considering  the  projection  of  x”  on  S(F).  | 

3.4  Remarks 

1 .  It  is  easy  to  see  that  Theorem  2.8  may  be  extended  to  the  nonlinear 
case  if  F  is  Lipschitzian.  In  this  case  ||M||  is  replaced  by  the  Lipschitz 
constant  of  F  in  (2.10). 

2.  Unfortunately,  from  a  computational  point  of  view  the  fixed  point 
methods  in  general,  and  those  considered  in  this  paper  in  particular,  are  not 
viable  methods.  They  are  extremely  slow  and  particularly  so  in  the  vicinity 
of  a  solution  point  since  the  step  sizes  taken  in  such  a  vicinity  are  extremely 
small.  Their  slowness  in  part  is  also  due  to  the  fact  that  they  do  not  utilize 
special  features  of  the  matrix  M  in  the  case  of  LCP(M,  g) .  Their  real 
utility  is  perhaps  in  generating  good  starting  points  for  fast  Newton-type 
algorithms.  However,  the  SOR  methods  are  much  faster  than  the  fixed  point 
methods  even  for  generation  of  starting  points. 
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